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For higher-derivative f{R) gravity where 7? is the Ricci scalar, a class of models is proposed 
which produce viable cosmology diflerent from the LambdaCDM one at recent times and satisfy 
cosmological. Solar system and laboratory tests. These models have both flat and de Sitter space- 
times as particular solutions in the absence of matter. Thus, a cosmological constant is zero in 
flat space-time, but appears effectively in a curved one for sufficiently large R. A 'smoking gun' 
for these models would be small discrepancy in values of the slope of the primordial perturbation 
- - - , power spectrum determined from galaxy surveys and CMB fluctuations. On the other hand, a new 

. problem for dark energy models based on f{R) gravity is pointed which is connected with possible 

3r ' overproduction of new massive scalar particles (scalarons) arising in this theory in the very early 

' Universe. 
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I. INTRODUCTION 



'■^ ' Continuing investigation of dark energy (DE) properties in the Universe (see the recent review [Jj for the definitions 
of what is usually called the effective DE energy density poE and pressure pde from the observational point of view) 
O has shown that its properties are very close to those of an exact cosmological constant A that has pA = —pA ~ 
^ A/SttG = const > 0} In particular, if wde = Pde/ Pde is assumed to be constant, then \wde + 1| < 0.1 (la 
^ • error bars) or even smaller, see Q for the analysis of the most recent observational data using different techniques.^ 
However, for more generic DE models with wde ^ const, the same analysis does not exclude varying wde including 
even temporal phantom behaviour of DE {wde < —1) at recent redshifts z < 0.3. The latter behaviour (in other 
words, breaking of the weak energy condition for DE), if confirmed by future, more exact data, may not be explained 
^ in the scope of physical DE models according to the terminology of [l| (e.g., the quintessence ones) and requires some 

kind of geometrical DE (otherwise dubbed modified gravity) , see also recent reviews Q . 
, Among geometrical DE models one of the most simplest ones is the f{R) class of gravity models with the Lagrangian 
psj ' density 



fiR) 



o: L = ^ + L„, (1) 

o 



16nG 

where Lm describes all non-gravitational kinds of matter including non-relativistic (cold) dark matter and the metric 
variation is assumed.'^ /(-R) rnay be an arbitrary function subjected, however, to some stability conditions discussed 
below. For f{R) = R — 2A, it reduces to the Einstein gravity with a cosmological constant. Thus, it contains the 
standard cosmological ACDM model as a particular case. However for f"(R) ^ 0, in addition to the massless spin-2 
' graviton, this class of models contains a scalar particle, dubbed scalaron in 5], which rest-mass is M^{R) — (3/"(i?))~^ 
in the WKB-regime \M'^\ R^, R^^^R'^'^ (here and below prime means differentiation with respect to an argument). 
It is neither a tachyon, not a ghost for f"{R) > in this regime. Finally, graviton is not a ghost if f'{R) > 0. All these 
properties can be easily obtained either directly, or (in the absence of Lm) using conformal equivalence of equations 
of motions for this class of models to those of the Einstein gravity interacting with a minimally coupled scalar field (p 
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^ The sign conventions are: the metric signature (H ), the curvature tensor R^^^i,^ = djyV^p — R^i, = so that the Ricci 

scalar R = R'^ > for the do Sitter space-time and the matter-dominated cosmological epoch; c = ?!, = 1 is assumed throughout the 
paper. 

^ Note that the assumption W£ie = const < and not equal to —1, —2/3, —1/3 is not very natural in the presence of non-relativistic 

matter since it requires DE models with rather specific potentials, see e.g. Q|. 
^ The Palatini variation of JTJ leads to completely different equations of motion, even the number of degrees of freedom (i.e. particle 

content) is not the same. Thus, two models with the same Lagrangian density |[T]| but different ways of variation should be studied and 

compared with data as two separate models; their similarity is illusory. 
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with some potential V^cj)) which form is uniquely determined by /(i?) in all points where f'{R) 7^ However, I 
shall not use this conformal equivalence below, in particular, because it may be misleading in the presence of other 
kinds of matter. 

In this paper Eq. ([T]) is considered as a purely phenomcnological effective Lagrangian density describing geometrical 
DE. However, in principle, it may arise either due to quantum fluctuations of all fields including gravity (as was 
supposed in early papers like [Hi), or as a result of reduction from higher dimensions to 4D in some variant of modern 
string/M-theory (see e.g. Q in this respect). 

First papers on cosmological models in f{R) gravity appeared already in 1969-1970 ^8,]. Then, among other results, 
this class of models with /(i?) = R + R^ /GM"^ plus some small non-local terms (which are crucial for reheating 
after inflation) was used to construct the first internally self-consistent cosmological model possessing a (quasi-)de 
Sitter (latter dubbed infiationary) stage in the early Universe with slow-roll decay, a graceful exit to the subsequent 
radiation-dominated Friedmann-Robertson- Walker (FRW) stage (through an intermediate matter-dominated one) 
and sufficiently effective reheating in the regime of narrow parametric resonance Q (see Q for more details). It is 
for this model that calculations of scalar (adiabatic) perturbations generated during inflation were first done [13] ■ 
Moreover, with all recent observational data taken into account, it still remains among viable cosmological models: 
as follows from the final results for both scalar and tensor perturbations [ll| (see also [l^, [3 [H)]), the model 
predictions for the slope of the primordial spectrum of scalar perturbations and the tensor/scalar ratio r are 
Hs-l = -2iV-i = -0.04 {N/50)-\ r = UN'"^ = 0.0048 {N/m)-'^ where N is the number of e-folds between 
the first Hubble radius crossing of the present inverse comoving scale 0.05 Mpc^^ and the end of inflation - well in 
agreement with the present data (of course, in the case of r we have an upper limit only). The only free parameter of 
this model - the scalaron rest-mass M - is determined from the normalization of the primordial scalar spectrum. If 
we take it from the best fit to the combined WMAP3-SDSS measurements [TH, then M = 2.8 x 10"^ (iV/50)"^Mpi 
where Mpi = l/\/G (by the way, this lies inside the range conjectured in [Tl|). 

Due to remarkable qualitative similarity between the present DE and primordial DE that supported inflation in 
the early Universe, all inflationary models may be applied to the description of the present DE, too, after changing 
numerical values of their microscopical parameters only.^ Actually, it had been done already for practically all models 
(sometimes in the inverse historical order). The same occurred to models based on f{R) gravity beginning from [l6| , 
and then it was proposed in [l3| to use f{R) models with f{R) diverging (or flnite but non-analytic) at i? for 
description of the present DE. However, after much agitation on this particular class of models, it was proven that 
they are either non-viable, or practically indistinguishable from the standard ACDM model, see [ij, [3| for such 
rather pessimistic conclusions as well as for extensive lists of publications on this topic. 

This does not completely close the way to construct a viable DE model in f{R) gravity observationally distinguish- 
able from the ACDM model, but suggests to abandon the hypothesis of divergence of f{R) aX R — as well at any 
other value of i?, and to return to the natural assumption that f{R) is regular in this point. Moreover, an interesting 
and intriguing possibility is /(O) = while / ^ R — 2A for i? ^ A. This behaviour corresponds to an effective 
cosmological constant existing in a sufficiently curved space-time but 'disappearing' in the flat one - that explains the 
title of the paper. In other words, in such a model the observed DE (close to the cosmological constant for sufficiently 
large R) is a purely curvature induced effect. It is totally unrelated to quantum vacuum energy in flat space-time 
that should be zero due to some other symmetry.^ However, the price to pay is that flat space-time becomes unstable 
with a characteristic time of the order of the present Universe age. 

Of course, it is much more difficult to construct a viable DE model in f{R) gravity as compared to general scalar- 
tensor gravity since the former contains only one arbitrary function while the latter has two functions and provides 
much place for viable DE models, see e.g. [20] for reconstruction of such models from different kinds of observational 
data and 21] for models which admit recent phantom behaviour of DE. Another source of problems is that it is rather 
non-trivial to satisfy laboratory and Solar system tests in /(i?) gravity since it formally represents the limiting case 
w = of scalar-tensor gravity where uj is the Brans-Dicke parameter [2^ , though sometimes this limit should be taken 
carefully (see [2^ for recent reconsideration). However, this may be considered even as an advantage of DE f{R) 
models since it makes easier to falsify them. So, in the next section a trial 3-paramctric form of f{R) is introduced 
which realizes the 'disappearing cosmological constant' possibility and behaviour of its FRW solutions is investigated. 
In Sec. 3 laboratory and Solar system tests, as well as dynamics of small perturbations are considered. Sec. 4 contains 
conclusions and discussion of problems and further tests of this model. 



This does not mean that the present and primordial DE should necessarily be the same kind of matter, like in the so called quintessential 
inflationary models. I am only speaking about possibility to use the same kind of theoretical models in both cases. 
^ So, this model easily realizes a possibility of 'degravitation of the cosmological constant' in the scalar sector which was recently proposed 
and tried to be achieved in the tensor sector in (Tj^l but faced with very serious technical problems due to non-locality of equations. 
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II. GEOMETRICAL DARK ENERGY MODEL AND BEHAVIOUR OF ITS FRW SOLUTIONS 

Field equations following from ([T]) can be written in the following Einsteinian form (though gravity itself is not the 
Einstein one): 

^m-^'^m^ = -8'^g(7;V) + t;(^^)) (2) 

where 

»^GT;^oe) ^ F'{R)R; - lF{R)d^^ + (V^V- - S;VpVP) F'{R) , EE f{R) - R (3) 

and follows from variation of Lm and satisfies the generalized conservation law T^.^^^^-) = separately (since 

the left-hand side of Eq. ^ and the right-hand side of Eq. 1^ satisfy this condition, too). There exists a subtlety in 
this representation that is discussed below. The trace of Eq. ([2]) reads 

3V^V^/' - Rf + 2/ = SttGT™ . (4) 

Constant curvature solutions (de Sitter ones for i? > 0) are roots of the algebraic equation Rf — 2/. 
Let us take f{R) in the following 3-parametric form: 

f{R)^R + XRjfl + ^\ -ij (5) 



with n,X > and Rq of the order of the presently observed effective cosmological constant. Then /(O) — (the 
cosmological constant 'disappears' in flat space-time) and R'j^ = is always a solution of Eq. ^ in the absence 
of matter, but /"(O) is negative - flat space-time is unstable. For \R\ ^ Rq, f{R) — R~ 2A(cx3) where the high- 
curvature value of the effective cosmological constant is A(oo) = XRq/2. The equation for de Sitter solutions having 
R = const = Ri = xiRq, xi > can be written in the form 

2((l + x2)"+l-l-(n+l)x?) ■ ^ ' 

Below, by xi I will mean the maximal root of Eq. ([5]). So, instead of specifying A, one may take any value of Xi 
and then determine the corresponding value of A. It follows from the structure of Eq. ^ that xi < 2A. Thus, the 
effective cosmological constant at the de Sitter solution A(i?i) = i?i/4 < A(oo). On the other hand, xi — > 2A in both 
limiting cases Xi fixed, n ^ 1 and Xi ^ 1, n fixed. In these cases the Universe evolution becomes indistinguishable 
from that in the ACDM model. 

Let us now consider the stability conditions 

/'(i?) >0 , /"(i?) >0 , R>Ri. (7) 

Note that they are imposed not in the whole space of solutions but only on a trajectory of the evolution of our Universe 
from very large and positive R in the past to i? = i?i in the infinite future. The quantum meaning of these conditions 
(graviton is not a ghost, scalaron is not a tachyon) has already been mentioned in Sec. 1. However, violation of 
these conditions in the course of purely classical evolution is undesirable, too. If f'{R) = for some finite R > Ri, 
a universe generically becomes strongly anisotropic and inhomogeneous at some finite moment of time (2^ . [25j (the 
same happens to the Einstein gravity -I- a non-minimally coupled scalar field, preventing Ge// from changing sign 
[26,]). In the point where f"{R) = 0, some weak singularity occurs which will be considered elsewhere. In terms of 
the conformal equivalence mentioned above, dR/d(j) diverges at this point. 

It can be shown that it is sufficient to satisfy the conditions (l7|) for i? = i?i and then they will be valid over the 
whole interval [Ri, 00). Correspondingly, this gives two necessary conditions for parameters of the form ([5|): 

(l-Hx?)"+i > l + (2n+l)a;2 , > l/(2n -f 1) . (8) 

To these inequalities, the condition of the stability of the future de Sitter stage has to be added. It follows from 
variation of Eq. ^ and reads (since the condition /"(i?i) > is already assumed to be satisfied) |i27i] : 



fiRi) > Rif'iRi) 



(9) 
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This condition is stronger than the first of inequaHties ([7]) at i? = so it substitutes it. For our model it produces 
the requirement: 

(1 + a;?)"+2 > 1 + (n + 2),t? + [n + l){2n + l)x\ . (10) 

It can be proven that it imphes the second of inequalities ([S]), too. So, it is sufficient to check this inequality only. In 
particular, for n = 1 it reads xi > -v/S that leads to A > 8/3\/3. In addition, the value of Xi saturating (|10p is also 
the point where X{xi) in Eq. ^ reaches minimum. 

Now turn to the model evolution at the matter dominated stage in the regime Rq. By construction, the model 
([5]) satisfies the conditions 

\F\^R, \F'{R)\^1, R\F"{R)\^1 (11) 

for i? ^ i?o. Due to this, it is possible to solve Eq. ^ iteratively in this regime and corrections to the standard 
Einsteinian (actually, even Newtonian) behaviour R = i?'^"-' = SirGTm (x ap pear to be small {a{t) is a FRW scale 
factor). Thus, the model does not possess the Dolgov-Kawasaki instability [2g|. These corrections are of two types - 
matter induced ones and free scalaron oscillations. The former ones follow from direct iteration of Eq. ([4|: 

R = + dR,„d + SRosc , <5i?md = (RF'iR) - 2F{R) - 3V^V^^^'(i?))^^^(„, (12) 

where is taken with respect to the unperturbed metric, too. For i? 3> -Ro, SRmd ~ const = — 2F(oo) = 2Ai?o = 
4A(oo) and ^''^GT'^^j-^^-^ « A(oo)(5J^. Thus, DE behaves as a positive cosmological constant at redsifts z ^ 1. 

Let us now return to the subtle point in the definition of the DE energy-momentum tensor ^ pointed above. It 
consists in the following: what is the constant G in it? For the model involved, G coincides with Ge//(oo) - the 
value of the Newtonian gravitational constant measured in a Cavendish-type experiment made in an environment 
with space-time curvature R ^ Rq. In the next section, we will discuss under what restriction on n this is achieved 
in laboratory experiments already made. Moreover, it coincides with Ge// at the radiation-dominated stage in the 
Universe during the period of the Big Bang nucleosynthesis (BBN) i ~ (1 — 100) s.^ So, BBN predictions remain 
unchanged in this DE model, too. Thus, such a choice of G in the left-hand side of ([3]) is very natural for this model. 
Taking another constant, say G{df / dR)^l^ where the subscript 'now' means the present moment, results in adding a 
'tracking' component to DE which is proportional to the Einstein tensor and may have any sign. That is why with 
our definition ([3]), there is no effect oi poE becoming negative at sufficiently large redshifts (with wde diverging at 
the moment when pjjE = 0) like that proposed in the recent paper [25|. 

More unexpected is the behaviour of the second small correction SRosc- It satisfies the equation 



i d f ^ d 
a^dt 



{a'^j^ (F"(i?(0)) SRosc)^ + SRosc = . (13) 



For R:^ Rq, the WKB-approximation may be used. Then the solution is 

with C = const. The integral here is just J M{R^^^)dt where M is the scalaron mass introduced in Sec. 1. At the 
matter-dominated regime: 

aoct2/3^ i?(o)^J_^ F" o^t^''+\ M(i?) oc t-2"-2 , (15) 

SRosc o^t-~^''~'^sm {const -t-^''-^) , ^ oc T sin (const • i^^n^^) . (16) 
Similar behaviour continues during the radiation-dominated stage: 

CX ^ M{R) CX t-3("+i)/2 , SRosc oc sin [const ■ t-(3«+i)/2) . (17) 



^ Note that at t ~ 100 s, after antimatter - positrons - annihilation, Tm is ~ 10 ^ g cm ^ only, similar to conditions on the Earth. 
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Thus, though oscillations of the scale factor remain small as t ^ 0, oscillations of R grow to the past and finally violate 
the assumption |5i?osc| *C Note that the energy of scalaron oscillations at that moment Posc ~ (SRosc)^ / GM"^ (R) 

is still much less than since R ^ for i? 3> i?o ■ Investigation of further FRW evolution of the model to the 
past is blocked by the stability problem: R can become less than i?i and even negative during oscillations, so the 
conditions ([7]) are violated. For rt > 1, it is possible to choose xi or A in such a way that f'{R) > for all real R. 
But f"{R) always becomes zero at i? = ±i?2 where R2 = i?o/\/2n- + 1 < -Ri and negative for \R\ < R2, and it is not 
possible to avoid this property without abandoning the assumption of the cosmological constant disappearance in flat 
space-time /(O) = 0. 

Therefore, first, to avoid M{R) becoming too large, say larger than Mp;, the function f{R) in (O has to be modified 
at R ^ 00. The simplest way is to add the term i?^/6M^ to ([5]) where the value of M (which will be the limiting 
value of M{R) for R 00) may be taken just that which is needed for the R + R^ infiationary model mentioned in 
Sec. 1. This term will be negligible for the present DE. 

Second, a new serious problem for DE models in f{R) gravity appears which has not been considered before: to 
avoid destroying radiation- and matter-dominated FRW stages, some mechanism in the early Universe should work 
to prohibit overproduction of scalarons. In the R + R^ infiationary model such mechanism does exist - gravitational 
creation of non-conformally-invariant particles and antiparticles (though not gravitons!) by oscillations of R. However, 
in the model ([5]) it is not possible to use it in full until the question how to evolve through the point where f"{R) — 
is solved. The only way to avoid this problem at all is to assume that |(5i?osc| < 87rGT„j just from the very beginning 
of early evolution of the Universe (or at least from the moment when the effective Lagrangian density (H]) becomes 
valid). In any case this makes the constant C in Eq. very small, practically zero at present and greatly reduces 
possible phase space for this DE model. 

Still, if this problem is solved and C may be put zero somehow, then a viable FRW background is obtained once 
the condition PH)) is satisfied. Moreover, since it follows from the present observational data that A(i?i) is not too 
much different from A(oo), one can choose xi close to 2A. Then the whole FRW evolution is analytically described 
by the first iteration formula (|12p where the background FRW metric is that of the ACDM model: 

aWocsinh^/^Qi/oi) , H^'^ = H, coil, (^-H,t^ , rW ^ ^^h^ + --^^^^ (18) 

where H = a/ a, H§ = A{oo)/3 — XRo/6. Knowing SRind, it is straightforward to obtain SH and 5a/ a. 

III. LABORATORY AND SOLAR SYSTEM TESTS AND DYNAMICS OF INHOMOGENEITIES 

The main problem with laboratory and Solar system tests of gravity for all f{R) DE models irrespective of the 
form of f{R) is that this class of models represents the limiting case w ^ of scalar-tensor gravity as was mentioned 
above. As a result, would scalaron be massless, an additional 'fifth force' would show itself in laboratory experiments 
and the first post-Newtonian parameters would have the values /3 = 1, 7 = 1/2 that is not admissible. However, this 
problem is qualitatively the same as it occurs for the string theory dilaton which in the massless low-energy limit 
corresponds to scalar-tensor gravity with u = —1 that is excluded, too. As is well known, that problem is solved by 
assuming that the dilaton is sufficiently massive. The same one has to assume about the scalaron in f{R) DE models. 
Namely, let us choose model parameters in such a way that for any laboratory or Solar system test of gravity having 
a characteristic scale L and made in an environment with some non-relativistic matter density pm, the scalaron mass 
satisfies the strong inequality M{R{pm))L ^ 1. Under this condition deviations from the equation R = SnGTm (and 
the Poisson equation, too) are small both inside and outside the Sun and other compact bodies. There have been 
hopes that this requirement may be circumvented using the so called "Chameleon" effect (s^l (see also [3l|). However, 
until recently it has not been shown to what extent can it be important in construction of viable /(i?) gravity models, 
see e.g. For the model ©, A/(p,„) oc p^^+i for i? > i?o with M(i?o) ~ iJ^^ - lO^^s cni-\ 

For the most recent and best Cavendish-type experiment [111, taking L « 50 /xm and p,„ « 10^^^ g cm^'^ (corre- 
sponding to a vacuum of ~ 10~^ torr achieved there), we obtain the sufficient condition n > 1 (rounded to a larger 
integer). In the case of light deflection or the Shapiro time delay by the Sun, the main contribution to 7 — 1 is from 
distances r ~ Rq « 7 x 10"'^° cm. If we take the Solar corona density at this distance (~ 10~^^ g cm~^) as p™, 
then already n > 0.5 would be sufficient. In the case 0.5 < n < 2 the condition M{R) r ^ 1 can be violated farther 
from the Sun. However, the scalar component of Solar gravitational field has been already screened by the Yukawa 
damping factor exp (— J M(r) dr) by that distance. Clearly more careful calculation is needed here with the detailed 
account of interplanetary matter profile in the Solar system. The same refers to such tests like Lunar laser ranging. 
In any case, all known laboratory and Solar system tests of gravity are certainly satisfied for n > 2 and probably this 
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condition may be softened up to n > 1. However, it will be shown below that there is no necessity in such softening 
due to limits from the large-scale structure of the Universe. 

Finally, let us turn to the evolution of weak inhomogeneities in the linear regime. It follows from equations for 
perturbations, either obtained for scalar-tensor gravity with the limit w — > taken in them, see e.g. [2u|, or from 
those directly derived for /(i?) theory ([l^, [H, [s^ and other papers), that the equation for density perturbations 
dm = Spm/ Pm in thc non-relativistic matter component (cold dark matter + baryons) during the matter-dominated 
stage has the form 

Sm + 2HSm - ^TrGeffPmSrn = (19) 

in the limit k ^ aH where k = |k| (spatial dependence exp(ikr) is assumed) and Ge// = G/f'{R) for k ^ M{R)a 
and Geff = AG/if'{R) for k > M{R)a {f'{R) « 1 for R > Rq). This increase of Ge// in 4/3 times is just how the 
'fifth force' due to additional scalar gravity shows itself at large scales and low matter density. 

Because of this, Sm grows as usually, oc t^/^, before the time moment tk when k = M{R)a, but after that this law 

changes to S,n oc that continues up to iA - the end of the matter-dominated stage when d ~ 0. Using (fT5|) . we 

obtain tk oc As a result, 5m{k) acquires an additional growth factor during the matter-dominated stage 

(compared to the ACDM model) proportional to 

✓ 33-1 2 
t A \ 6 3 ^33-5 

— CX /fc^T3^ . (20) 

tk J 

This additional increase occurs at small redshifts and is not seen in CMB fluctuations (apart from some features in 
a few low multipoles which are bounded by cosmic variance). E.g., for n = 2, k/anowHnow — 300 and z\ — 0.7, 
the redshift Zk{tk) = 2.5. As a result, there arises some discrepancy between values of the slope Us of the primordial 
power spectrum determined from galaxy surveys on one side (assuming the standard evolution of perturbations) and 
CMB fluctuations on the other: 



2(3n-h2) ^ ' 

For comparison, An, is equal to 0.074 for n = 1 and 0.047 for n = 2. Note that the limit n corresponds not to 
the ACDM model but to F{R) cx Ini? at large R. 

At present no such discrepancy is seen, see e.g. so we may conservatively bound Ans < 0.05 that leads to 

n > 2 in ([5]). Of course, a more exact numerical calculation is needed for larger values of n since then the values of 
Zfe lie rather close to unity, so the formula (|2ip becomes too approximate. 



IV. CONCLUSIONS AND DISCUSSION 



Thus, in contrast to numerous unsuccessful previous attempts to construct a viable DE model in ,f{R) gravity using 
a function / divergent or non-analytic at i? = 0, it appears possible to achieve this goal with a regular f{R) satisfying 
the condition /(O) = which means the absence of a 'bare' cosmological constant in flat space-time. On the other 
hand, DE in this model behaves itself as an effective cosmological constant for large R if F{R) = f — R ^ const at 
i? — > oo. This model passes laboratory and Solar system tests of gravity if its parameter n is suflficiently large (n > 2 
seems to be the sufficient, but probably not necessary condition), though analysis of gravitational radiation from 
double pulsars may add some new restriction. But it is clear already that limits from large-scale structure arising 
due to the anomalous growth of linear perturbations at recent redshifts are more critical and lead to a stronger limit 
on n. Just the opposite, any discrepancy An, between values of the slope of the primordial perturbation spectrum 
obtained from galaxy and CMB data may serve as a strong argument for such model. Then Eq. pT|) directly relates 
Aug to n. 

However, more deep theoretical analysis of this class of models has uncovered its new serious problem not considered 
before: how to avoid an overabundance of new scalar particles arising in f{R) gravity (dubbed scalarons) which can 
be generated in the very early Universe. Mathematically this means that the coefficient C in Eq. (fT4|) should be 
practically zero at present, i.e. almost the whole degree of freedom has to be suppressed. Otherwise, a FRW solution 
in this model cannot have sufficiently long radiation- and matter-dominated stages because it hits the weak singularity 
where f"[R) — 0. Note that this difficulty is of a more subtle type than those of hnear stability considered before. 
It is rather a problem of a measure of initial conditions in the early Universe leading to the standard cosmological 
evolution almost up to the present time. Anyway it remains an important topic for further study. 

While this paper was being prepared for publication, two papers (s^ . Issj appeared where similar DE models in 
/(i?) gravity possessing the property /(O) = were proposed. Our model is closer to that in |34l] . 
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